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ABSTRACT 


This  report  reviews  a generalized  formulation 
of  the  steady-state  boundary-value  problem  for 
scattering  of  infinite  plane  waves  by  an  arbitrary 
closed  rigid  surface  immersed  in  an  infinite  fluid. 
The  normal  velocity  distribution  generated  over  the 
closed  surface  by  scattering  of  plane  waves  can  also 
be  interpreted  as  the  boundary  condition  of  an  equi- 
valent steady-state  radiation  problem.  The  numerical 
solution  of  rigid  surface  scattering  problems  is 
therefore  obtainable  by  a simple  extension  of 
capabilities  of  the  XWAVE  program.  The  additional 
data  required  by  XWAVE  for  rigid-surface  scattering 
applications  and  several  sample  calculations  are 
presented. 


FORMULATION  OF  BOUNDARY -VALUE  PROBLEM  FOR  RIGID-BODY  SCATTERING 

Part  I^  of  this  documentation  described  a numerical  method  for 

obtaining  the  radiated  pressure  field  external  to  the  surface  of  a 

structure  vibrating  in  an  ideal,  infinite  fluid.  This  problem  is 

mathematically  posed  by  the  wave  equation  with  attendant  boundary  conditions 

1 2 3 

at  infinity  and  on  the  structural  surface.  It  is  well  known  ’ that 
boundary-value  problems  such  as  the  one  describing  sound  pressure  from 
submerged  vibrating  surfaces  can  also  be  interpreted  as  formulating  the 


Henderson,  F.M. , "A  Guide  to  Use  of  the  XWAVE  Program:  Part  I - Radiated 
Pressures  From  Vibrating  Structures,"  DTNSRDC  Report  77-0041  (Jun  1977). 

2 

Chertock,  G. , "Integral  Equation  Methods  in  Sound  Radiation  and 
Scattering  from  Arbitrary  Surfaces,"  NSRDC  Report  3538  (Jun  1971). 

3 

Junger,  M.C.,  and  D.  Feit,  "Sound,  Structures  and  Their  Interaction," 

MIT  Press,  Cambridge,  Massachusetts,  and  London,  England  (1972). 
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sound  field  produced  when  trains  of  pressure  waves  impinge  on  and  are 
scattered  by  structural  surfaces. 

To  obtain  this  alternative  formulation,  the  vibrating  structure 
considered  in  the  earlier  report^  is  replaced  by  a rigid-body  surface  (see 
Figure  1)  on  which  a time-independent  succession  of  infinite  plane  waves 
impinges. 


Figure  1 - Rigid-Body  Scatterer  Immersed  in  an  Infinite  Fluid 

It  is  of  practical  interest  to  consider  plane  wave  incidence  since  in 
applications  such  as  active  sonar  detection  problems  the  wave  source  is 
sufficiently  far  from  the  scattering  obstacle  that  the  radius  of  curva- 
ture of  the  wave  fronts  near  the  obstacle  can  be  considered  essentially 
infinite.  In  addition,  the  decrease  in  pressure  from  geometric  spreading 
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of  the  incident  wave  over  the  scattering  surface  can  be  neglected. 

The  total  pressure  (above  ambient)  in  the  field  is  the  sum  of  the 
incident  and  scattered  wave  pressures, 


p(z)  = p.(z)  + Pg^(z) 


(1) 


where  z_  denotes  field  points  and  p^^,  in  the  notation  of  Junger  and  Feit,' 


denotes  field  pressures  scattered  from  rigid  (infinite  impedance) 
boundaries.  Since  the  total  pressure  p must  satisfy  the  wave  equation, 
the  scattered  pressure  must  also  satisfy  it. 


v2p  (z)  + k^p  (z)  = 0 


(2) 


k is  the  wave  number  of  the  incident  pressure  wave;  k = w/c  in  which 
u)  is  the  angular  frequency  of  the  incident  pressure  wave 
c is  the  speed  of  sound  in  the  fluid 

The  presence  of  the  rigid  boundary  requires  the  sum  of  the  fluid  particle 
velocities  in  the  incident  and  scattered  waves  at  any  point  on  the  body 
surface  to  be  zero. 


sp^ls) 


an 


3Pc„(s) 

iu)pv.(s)  = iojpv^  (s)=- 
1 s°° 


8n 


(3) 


where,  referring  to  Figure  1, 

s is  a point  on  the  closed  surface  S of  arbitrary  shape 
n is  the  direction  normal  to  the  structural  surface  at  s 
v^(s)  is  the  incident  fluid  particle  velocity  normal  to  the 
structural  surface  at  s 


v^_^(s)  is  the  scattered  fluid  particle  velocity  normal  to  the 


structural  surface  at  s 


p^(s)  is  the  incident  pressure  on  the  structural  surface  at  s 


p^  (s)  is  the  scattered  pressure  on  the  structural  surface  at  s 

Soo 


p is  the  fluid  density 
The  boundary  condition  at  infinity  is 


.iklzl 


(4) 
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I 
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t 
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where  ^ denotes  a point  in  the  fluid. 

Equations  (2),  (3),  and  (4)  then  give  the  alternative  boundary- value 
problem  to  be  solved  for  rigid  body  scattering.  When  these  equations 
are  compared  with  their  counterparts  for  a vibrating  body,  (Equations  (1), 
(2),  (3)  of  DTNSRDC  Report  77-0041)  it  is  seen  that  the  only  distinction 
between  the  radiation  problem  and  the  rigid-body  scattering  problem  is 
the  manner  in  which  the  surface  boundary  condition  is  interpreted. 


GENERAL  FORMULATION  FOR  THE  BODY  SURFACE  BOUNDARY  CONDITION 


In  order  to  perform  rigid-body  scattering  calculations  using  XWAVE, 
a suitable  representation  for  the  normal  surface  boundary  condition 
Vsoo(s)  arising  from  an  incident  plane  wave  must  first  be  obtained.  Since 
a general  direction  cf  incidence  is  to  be  specified,  it  is  convenient  to 

3 

utilize  vector  wave  number  notation.  With  this  notation  the  steady-state 
incident  pressure  field  is  given  by 

ik.(-)2 

p.(z)  = P.e  (5) 


where  is  defined  as  a vector  with  magnitude  equal  to  and  having  x-, 
y-,  z-components  corresponding  to  the  direction  of  wave  incidence,  and  z 
is  a vector  denoting  a field  point.  Taking  the  partial  derivative  of 
p^il)  in  the  direction  normal  (outward)  to  the  surface  S at  s (z^  = s) 
gives 


8Pi(s) 

3n 


9P  .e 
1 


ik,(-)s 


-1 


an 


= ^-vP^.e 


ik.(-)s 


-1 


= i ri-k.^P.j  e 


ik,(-)s 


-1 


(6) 


where  ^ denotes  a unit  vector  normal  to  S at  s and  pointing  to  the 
exterior  of  S.  Substituting  this  result  into  Equation  (3)  then  yields 

ik.  (Os 

in.  k.P.e  = iu)pv^^(s)  (7) 

Since  the  XWAVE  formulation  uses  the  nondimensional  forms  of  pressure  and 


i 


4 


velocity 


I 


P = P/pCVq 
7 = v/Vq 


(8) 


where  is  an  arbitrary  velocity,  dividing  both  sides  of  Equation  (7)  by 
ipcvQ  gives 

_ iki{-)s 

n-kP.  e =k.v(s)  (9) 

— -1  1 1 S” 


and  finally  the  desired  expression  for  the  surface  normal  velocity 
resulting  from  time-harmonic  infinite  plane  waves  incident  on  S from  a 
general  direction. 


V,  (s)  = 

S°o ' ' 


J1  ■ 


ik,(-)s 


k. 

1 


(10) 


Since  XWAVE  is  based  on  a numerical  solution  to  the  wave  equation,  a 
finite  representation  for  v^^(s)  is  used,  as  in  the  case  of  the  radiation 
calculations,  in  which  its  values  are  defined  only  for  points  of  a 
discretized  model  of  the  scattering  surface  S (see  Figure  2 of  DTNSRDC 
Report  77-0041). 


XWAVE  DATA  FOR  THE  RIGID-BODY  SCATTERING  PROBLEM 

Equation  (10)  shows  that  the  data  required  to  calculate  v^_^  are: 

1.  X-,  y-,  z-coordinates  of  ^ 

2.  X-,  y-,  z-coordinates  of  vector  wave  number 

3.  X-,  y-,  z-coordinates  of  surface  point  s 

4.  magnitude,  P.,  of  nondimensional  incident  pressure 

5.  magnitude,  k^ , of  1^^ 

Since  items  1)  and  3)  are  already  included  as  part  of  XWAVE's  "surface 
geometric  data",^  additional  input  facility  is  needed  only  for  the 
remaining  items.  To  accommodate  items  2)  and  4)  a new  data  card  is 
introduced  with  format  as  follows: 
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• "INCIDENT  PLANE  WAVE  DATA"  CARD 


I 


Col umns 
1-8 


11-18 

21-28 

33-40 


Contents 

k. 

-1 

X 

Description 
x-component  of  1^^ 

k. 

1 

y-cc  lonent  of  k. 

1 

1 

y 

k. 

-1 

z-component  of  1^^ 

z 

P. 

1 

Magnitude  of  nondimensional  incident 

pressure  wave,  P^.  = p^/pcv^ 


Item  5)  is  entered  in  the  space  allotted  to  k (for  radiation 
applications)  on  the  "Miscellaneous  Data"  CardJ 

DATA  INPUT  FORM  (1)  (see  Part  I of  this  documentation)  is  revised  as 
shown  in  Appendix  A to  include  the  incident  plane  wave  data  card  type  and 
to  indicate  its  position  in  the  XWAVE  data  deck. 

Two  of  XWAVE' s program  options  (see  Part  I,  section  entitled  DATA 
FORMATS)  are  augmented  as  follows: 


• "PROGRAM  OPTIONS"  CARD^ 


Col umns 
13-16 


17-20 


Contents 


0005 


0000 

0001 


Description 

0PT4:  Selector  for  surface  velocity 
di stri buti on 

Velocity  distribution  is: 

Generated  for  plane  wave  incidence 

upon  the  body  surface.  (Scattering 
applications  only) 

0PT5:  Selector  for  radiation  or 
scattering  applications 
Radiation  calculations 
Scattering  calculations 


The  other  options  and  types  of  data  as  well  as  data  configurations 
to  be  used  for  scattering  applications  are  the  same  as  those  previously 
described  for  use  with  radiation  calculations  (see  DATA  FORMATS^). 
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SAMPLE  V.ALCULATIONS 

CALCULATION  OF  FAR-FIELD  SCATTERED  PRESSURES  FROM  A RIGP  CYLINDER 
ILLUMINATED  BY  A PLANE  WAVE  NORMAL  TO  THE  CYLINDER  AXIS 

3 

An  approximating  analytic  result  is  used  for  comparison  with  the 
numerical  solution  for  this  problem.  The  analytic  approach  extends  results 
obtained  for  an  infinite  rigid  cylinder  that  scatters  a plane  wave 
incident  from  the  ((>=11  direction  (spherical  coordinates)  by  restricting  the 
scattering  acceleration  boundary  condition  of  the  infinite  cylinder  to  a 
finite  length  2L,  Figure  2,  and  by  taking  the  acceleration  distribution 
equal  to  zero  for  |z|  > L. 


0 » 0 


0 = 


Figure  2 - Approximate  Representation  for  Finite  Cylinder  Scattering  of 
Plane  Wave  Incident  Normal  to  Cylinder  Axis 


Ji-r 

A 


The  acceleration  for  IzUL  is' 


w(z,<t>)  = 


kP. 

1 


Z e„i”j'(ka)  cos  ncf. 
n=0  " " 


(11) 


with  the  following  notation'^: 


e is  the  Neumann  function,  =1  for  n=0;  2 for  n>0 
n 


J^(x)  is  the  cylindrical  Bessel  function  of  the  first  kind 


a is  the  cylinder  radius 


The  far-field  scattered  pressure  field  is  expressed  as' 


1 W R I 

2Le  P.j  (kL  cos  0)  e J (ka)cos  n4) 

Ps»(^’®’^)  = TT  R sin  0 H'  (ka  sin  0) 

n=0  n 


(12) 


with  the  following  notation  ; 

R,0,(|>  are  the  spherical  coordinates  of  field  point 
sin  X 


jo(x)  = 


H (x)  is  the  cylindrical  Hankel  function  of  the  first  kind 
n 


Before  Equation  (12)  is  evaluated,  it  is  convenient  to  modify  it 
to  make  the  results  compatible  with  XWAVE  calculations.  Transforming  to 
nondimensional  pressures  by  dividing  through  by  pcv^  and  using  Equation 
(8),  gives 


2Le^*^^P.  j„(kL  cos  e)  e j'(ka)cosn$ 
" TT  R sin  0 H|I|(ka  sin  0) 


(13) 


The  form  of  the  XWAVE  solution  for  far-field  pressure  (Equation  (10)  of 
DTNSRDC  Report  77-0041)  is 


. 4tt  I z I 

f'i)  ^F?<i>THTT 

e w 


(14) 


In  this  expression,  p(z)  denotes  nondimensional  pressure  at  z,  where  z is 


the  far-field  point  position  vector.  For  the  case  of  scattering. 


Ps_^(^)  will  correspond  to  p^^(R,9,(t>)  in  Equation  (13).  The  following 
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form  of  Equation  (13)  is  thus  compatible  with  XWAVE : 


F (z)  = ^ p (z) 
s®  k s®  V 


k 


(R,0,<l>) 


4tiR 

ikR 

e 


8iLP.  jp,(kL  cos  e)  c j'(ka)cos  n4> 

= - k s"n  e H''  ka  ?i-n  * fs-'*-*'  >’5> 

n=0  n 

To  illustrate  the  evaluation  of  Equation  (15)  the  following  data 
were  used:  cylinder  radius  a=l ; cylinder  half-length  L=2;  k=l ; normalized 
incident  pressure  P^=l- 

A FORTRAN  program  was  written  to  perform  the  far-field  pressure 

calculation.  This  program  utilized  a subroutine  COMBES  (PS-582A)  from  the 

Boeing  Math  Science  Library  to  compute  the  cylindrical  Bessel  functions 

of  first  and  second  kinds,  J (x)  and  Y (x),  for  this  evaluation.  These  in 

4 n n 

turn  were  used  to  obtain  : 

J'(ka)  = -J  ,,(ka)  + J (ka) 
n ' n+1 ' ' ka  n'  ' 


H^(ka  sin  e)  = J^(ka  sin  e)  + i (ka  sin  o)  (16) 

H^(kasin0)  = (ka  sin  0)+i  'n.l  (ka  sin  e)]  + ^ 

+ i Y (ka  sin  0)1 

Far-field  pressures  were  computed  around  a polar  great  circle  path  (in  a 
plane  perpendicular  to  the  incident  wave)  and  an  equatorial  great  circle 
path  (in  a plane  parallel  to  the  incident  wave)  on  the  surface  of  a large 
sphere  (R  ->  ®)  centered  about  the  cylinder.  The  results  are  given  in 
Table  1. 


Hildebrand,  F.B.,  "Advanced  Calculus  for  Appl ications ,"  (Prentice-Hall 
Inc.,  Englewood  Cliffs,  New  Jersey,  1962).  (See  in  particular  Chapter  4, 
Equations  (92),  (111).) 
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TABLE  1 - FAR-FIELD  SCATTERED  PRESSURES  FROM  FINITE.  RIGID  CYLINDER 
WITH  INCIDENT  PLANE  WAVE  NORMAL  TO  CYLINDER  AXIS  (ANALYTIC  SOLUTION) 


1 1/  1 ^5(9  ><t>) 

* 

max 

9° 

Polar  Circle  ({1° 

9° 

Equatorial  Circle 

0 

0 

.3397  90 

0 

.3397 

5 

.2932 

5 

.3389 

10 

.2436 

10 

.3400 

15 

.1922 

15 

.3451 

20 

.1432 

20 

.3548 

25 

.1102 

25 

.3694 

30 

.1177 

30 

.3887 

35 

.1642 

35 

.4126 

40 

.2274 

40 

.4404 

45 

.2956 

45 

.4713 

50 

.3639 

50 

.5044 

55 

.4299 

55 

.5383 

60 

.4920 

60 

.5717 

65 

.5484 

65 

.6031 

70 

.5977 

70 

.6312 

75 

.6382 

75 

.6545 

80 

.6683 

80 

.6720 

85 

.6870 

85 

.6829 

90 

.6933 

90 

.6866 

95  \ 

95  ) 

; 

Symmetry  (85°-0°) 

Symmetry  (85°-0°) 

180  ' 

180  ^ 

180 

180 

.3397  270 

180 

175 

.3821 

• 

Symmetry  (180°-0°) 

170 

.4265 

0 

165 

.4756 

160 

.5295 

155 

.5865 

150 

.6442 

145 

.7001 

140 

.7523 

135 

.7998 

130 

.8421 

125 

.8793 

120 

.9112 

115 

.9383 

no 

.9604 

105 

.9776 

100 

.9900 

95 

.9975 

■Ar 

90 

1 .000 

iFc 

(0>4>)L.v  = ■'^•8013 

85  \ 

' s 

'max 

I 

Symmetry  {95°-180°) 

0^ 

t 


In  order  to  obtain  the  numerical  solution  of  this  problem,  the 
finite  cylinder  of  length  2L  (see  Figure  2)  is  referenced  to  a Cartesian 
frame  as  shown  in  Figure  3, 


X 


Figure  3 - Plane  Wave  Incidence  Normal  to  Axis  of  Finite  Cylinder 

and  the  incident  wave  front  (Figure  2)  is  reoriented  to  approach  from 
the  y'*’  direction  (ij)  = 90°).  The  latter  revision  transforms  the  xz-plane 
symmetry  of  the  surface  velocity  boundary  condition  resulting  from  x" 

(or  (|i  = it)  incidence  to  a particular  pattern  of  yz-plane  symmetry 
compatible  with  XWAVE's  data  input  generator. 

As  a result  of  yz-  and  xz-plane  symmetry  in  the  cylinder  geometry 
and  velocity  boundary  condition,  specification  of  the  acoustic  model  for 
only  one  quarter  of  the  cylinder  surface  is  sufficient  to  enable  the  data 
generator  to  establish  the  total  surface  model.  Figure  4 indicates 
subdi visioning  of  the  quarter  surface  into  "regions"^  (numbered)  for 
subsequent  generation  of  surface  acoustic  elements. 

The  data  which  specify  the  surface  regions  and  the  distribution  of 
elements  over  each  region  (and  the  total  body  surface  by  reflection)  are 
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No,  Basic  Elements 
z:  n X m 


No.  Effective  Elements 
Over  Entire  Surface 
4(z  n X m 


MODEL  A 


MODEL  B 


162 

648 

270 

1080 

* Ref.  Figure  4 

Numerical  results  for  the  surface  modelings  given  in  Table  2 are 
summarized  in  Table  3 and  the  results  for  the  second  model  are  compared  in 
Table  4 and  Figure  5 with  results  for  the  approximating  analytic  solution 
(Table  1).  The  data  for  XWAVE,  prepared  on  the  program's  input  forms, 
and  computer  output  for  the  second  calculation  (Model  B)  are  presented  in 
Appendix  B. 
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IMDLL  o - rMK-riLLU  DLMMtr\LU  rr\Loour\Lo  rrwjn  riiiiiL,  iMUiu 

WITH  PLANE  WAVE  NORMAL  TO  CYLINDER  AXIS  (XWAVE  SOLUTION) 


W 
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TABLE  4 - PERCENTAGE  DIFFERENCES  BETWEEN  NUMERICAL  AND 
ANALYTIC  SOLUTIONS  FOR  FIRST  SAMPLE  PROBLEM 


Polar  Circle 


Equatorial  Circle 


0° 

% Differences* 

aO 

0° 

% Differences* 

0 

10.9 

90 

0 

10.9 

5 

15.5 

5 

11.7 

10 

24.0 

10 

12.9 

15 

39.3 

15 

13.7 

20 

65.4 

20 

13.9 

25 

93.1 

25 

13.4 

30 

70.9 

30 

12.2 

35 

26  J 

35 

10.3 

40 

35.6 

40 

8.06 

45 

5.14 

45 

5.54 

50 

7.36 

50 

2.91 

55 

6.84 

55 

0.37 

60 

5.30 

60 

-1.96 

65 

3.50 

65 

-4.01 

70 

1.81 

70 

-5.75 

75 

0.41 

75 

-7.12 

80 

0.64 

80 

-8.10 

85 

1 .27 

85 

-8.68 

90 

1.47 

90 

-8.88 

10.9 

8.98 

7.39 

5.47 

3.32 

1 .33 
-0.28 
-1.31 
-1 .82 
-1.90 
-1.70 
-1.34 
-0.94 
-0.59 
-0.29 
-0.072 
-0.071 

0.14 

0.20 


* Numerical  Result-Analytic  Result 
= A^^^y■t^c  RTsuI t " 100 
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.•jA'  . 


Figure  5 - Far-Field  Scattered  Pressures  from  Finite,  Rigid  Cylinder  with 
Incident  Plane  Wave  Normal  to  Cylinder  Axis.  Comparison  of 
Analytic  and  Numerical  Solutions 


Figure  5a  - Pressures  Along  the  Polar  Circular  Path  on  a Spherical 

Surface  at  Infinity 


Figure  5 (Continued) 


A 


I 


I 


-/7e 


e « 90° 


XWAVE  SOLUTION 


Figure  5b  - Pressures  Along  the  Equatorial  Circular  Path  on  a 
Spherical  Surface  at  Infinity 
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CALCULATION  OF  FAR-FIELD  SCATTERED  PRESSURES  FROM  A RIGID  CYLINDER 
ILLUMINATED  BY  A PLANE  WAVE  PARALLEL  TO  THE  CYLINDER  AXIS 


The  second  sample  problem  illustrates  the  facility  with  which  the 
direction  of  the  incident  wave  and  hence  the  velocity  boundary  condition 
for  scattering  can  be  varied  by  using  the  wave  vector  concept  as 
implemented  in  XWAVE. 

Instead  of  a plane  wave  incident  normal  to  the  axis  of  a cylinder 
(Figure  3),  a wave  propagating  parallel  to  the  z-axis  from  the  z direction 
is  used,  as  shown  in  Figure  6.  If  the  wave  number  is  assumed  to  be  2 for 


this  case,  the  components  of  the  wave  number  vector  are  k.  =0,  k^. 

=2.  Although  any  direction  in  space  could  be  as  readify  spec 

fr5m  the  general  expression  for  components,  ii^=2a,  M =2s,  ]ii,=2 

A y z 


=0, 

■^f  ied 
Y , where 


a,  B,  Y are  direction  cosines  of  k ^ , the  wave  directed  parallel  to  the 


z-axis  generates  a radially  symmetric  surface  velocity  boundary  condition 


which  in  turn  leads  to  shorter  computer  running  times  for  the  problem. 


The  portion  of  the  cylinder  surface  specified  to  XWAVE's  data 
generator  and  the  subdivision  of  this  surface  into  regions  used  by  the 
program  in  establishing  surface  element  models  are  sketched  in  Figure  7. 


X 


Figure  7 - Designation  of  Surface  Regions  on  Portion  of  Finite 
Cylinder  for  Second  Sample  Problem 

Table  5 summarizes  three  surface  models,  the  initial  one  and  two 
refinements,  for  which  calculations  were  made.  In  all  cases  the  far-field 
pressures  were  calculated  at  2-degree  intervals  to  insure  reasonably  good 
resolution  of  the  lobes  which  appear  in  the  pressure  profile.  Data 
preparation  and  computer  output  for  the  initial  calculation  (Model  1)  are 
given  in  Appendix  B. 

Table  6 compares  results  obtained  for  the  three  models.  The  initial 
surface  model  results,  which  differ  by  at  most  2.11%  from  the  second  model 
results  and  2.94%  from  the  third  model  results,  indicate  that  the  initial 
model  was  already  a fairly  good  representation  of  the  cylinder  surface. 

Figure  8 sketches  the  scattered  pressure  profile  predicted  by  th ? 
third  model . 
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V 

A 


TABLE  5 - SURFACE  MODELS  OF  FINITE  CYLINDER 
FOR  SECOND  SAMPLE  PROBLEM 


MODEL  3 


TABLE  6 - PAR-FIELD  SCATTERED  PRESSURES  FROM  FINITE,  RIGID  CYLINDER 
WITH  INCIDENT  PLANE  WAVE  PARALLEL  TO  CYLINDER  AXIS.  COMPARISON 
OF  XWAVE  RESULTS  FOR  THREE  SURFACE  MODELS. 


IF  (9)1 

e 

Model  1 

Model  2 

Model  3 

e 

Model  1 

Model  2 

Model  3 

0 

8.07 

8.21 

7.91 

92 

4.01 

4.04 

4.00 

4 

7.99 

8.13 

7.83 

96 

4.55 

4.58 

4.54 

8 

7.76 

7.89 

7.59 

100 

4.66 

4.68 

4.66 

12 

7.37 

7.50 

7.21 

104 

4.34 

4.35 

4.35 

16 

6.85 

6.96 

6.69 

108 

3.64 

3.64 

3.66 

20 

6.22 

6.31 

6.05 

112 

2.71 

2.68 

2.72 

24 

5.52 

5.60 

5.37 

116 

1 .84 

1 .79 

1 .83 

28 

4.86 

4.91 

4.73 

120 

1.79 

1.78 

1.74 

32 

4.37 

4.39 

4.26 

124 

2.72 

2.75 

2.64 

36 

4.20 

4.19 

4.12 

128 

3.89 

3.95 

3.81 

40 

4.38 

4.37 

4.34 

132 

4.98 

5.05 

4.90 

44 

4.81 

4.80 

4.79 

136 

5.85 

5.93 

5.78 

1 48 

5.29 

5.30 

5.28 

140 

6.47 

6.55 

6.41 

52 

5.65 

5.68 

5.64 

144 

6.82 

6.90 

6.78 

56 

5.75 

5.80 

5.74 

148 

6.94 

7.00 

6.91 

60 

5.53 

5.58 

5.51 

152 

6.87 

6.92 

6.86  i 

64 

4.94 

5.00 

4.92 

156 

6.66 

6.69 

6.66 

1 68 

4.02 

4.07 

4.00 

160 

6.39 

6.40 

6.39 

72 

2.84 

2.87 

2.82 

164 

6.09 

6.10 

6.10 

76 

1.56 

1.58 

1 .55 

168 

5.83 

5.82 

5.84  : 

o 

CO 

0.98 

0.960 

0.967 

172 

5.62 

5.61 

5.63  1 

84 

1.94 

1.94 

1 .94 

176 

5.50 

5.48 

5.50 

88 

3.10 

3.11 

3.09 

180 

5.45 

5.44 

5.45  ' 
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SCATTERING  APPLICATIONS 


title  DATE  

PROBLEM  SHEET  OF 


XWAVE  DATA 


PROBLEM  TITLE  

13 

CASE  TITLE  

13 


DIMENSIONS  FOR  XWAVE 


01 

02 

03 

04 

Ob 

06 

07 

08 

09 

DIO 

011 

012 

013 

014 

— 

1 

t 

5 

9 

13 

1 7 

31 

2b 

29 

33 

37 

41 

45 

49 

53  5( 

MISCELLANEOUS  DATA 

k. 

1 

( 

u 

H iREALl 

H (IMAG) 

ITER 

ATlON 

LIMIT 

CONVERGENCE 

CRITERION 

’ 

9 

17 

2b 

33 

41 

49 

53 

60 

PROGRAM  OPTION^; 

0PT3 

OPT4 

OPTb 

0PT6 

OPT? 

OPTS 

OPT9 

OPTIO 

1 5 9 13 W ?1 ^ ~79 T5 77  AO 

INCIDENT  PLANE  WAVE  DATA  


aaI 

0 

c>_ 

1  

Ia/w\ 

[ 

1 8 

1 11  18  ?1  ?8 

i 57  JT! 

SYMMETRY  OPTIONS  

PLANE  Pi  ANE  LONG,  radial 

"i  6 9 T3 


DATA  INPUT  FORM  (1) 


i 
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APPENDIX  B 

INPUT  DATA  AND  COMPUTER  OUTPUT  FOR 
SAMPLE  CALCULATIONS 


SCATTERING  APPLICATIONS 


title sRnpke ppoBkEn  / 

PROBLEM  


DATE  

SHEET  / OF  3 


XWAVE  DATA 


PROBLEM  TITLE 


rf  cykfNDEJ{^  ^ />  length*^ 


13 


CASE  TITLE 


A/w^AA/^/w^A^^/{i&/0"^ooy  scfirr£/UN&  of  firhe  mit  F^nyjf)  oif^icTion 


13 


D1 

D2 

D3 

D4 

Ob 

rns 

OOiO 

oVo 

aoo;^ 

60/S 

00  Jl 

ooof 

’ 

b 

9 

13 

17 

?1 

25 

29 

33 

37 

41 

45 

49 

53  56 

MISCELLANEOUS  DATA 

k. 


ITER 

w HMAPj  ation  convergence 

H IIMAGI  CRITERION 


1 

PROGRAM  0 

OPl  1 OPT2 

9 

17 

25 

33 

41  49  53  60 

OPT3  0PT4  0PT5  OPT6  OPT7  OPTS  OPT9  OPTIO 

oooooooi 

Ooo)  OOOl 

O 

o 

OOOO  OOOl 

1 5 9 13 

INCIDENT  PLANE  WAVE  D 
k,  k. 

■y 

17  21  25  29  33  37 

<^TA 

-i  p • 

i 1 

0.0 

-J.o 

0.0 

OC/)Ci 

1.0 

t 0 

11  18  21  20 

33  TC 

SYMMETRY  OPTIONS  

'ii  %i 

PLANE  plane  long.  RADIAL 


OOCZ\0O0O  0001  ooool 


TITLE 

PROBLEM 


PIECEWISE  CONICAL  SHELL  SURFACE 


f\OBLEM 


DATE  

SHEET  X.  OF . 


XWAVE  DATA 


NUMBER  OF  REGIONS  . 


\m 


1 4 

REGION  EXTENT  AND  MODELING  

r, 

02 

^2 

n 

m 

SIGN 

190.0 


9 17 

normal  surface  velocity 

V.  (REAL!  V_  (IMAG.I 


49  53  57 


1^91 


1 8 11  18 
REGION.  EXTENT  AND  MODELING 

'1  ^7  ^1 

»2 

^2 

n 

m 

SIGN 

1.0  1.0 


1 9 17 

NORMAL  SURFACE  VELOCITY 
V„  (REAL)  V„  (IMAG.) 


/?ao 


\m\ 


OOI^ 


49  53  57 


IWI 


1 8 11  18 

REGION  EXTENT  AND  MODELING  - 

M ^7  ^1 

«2 

^1 

*2 

n 

fT> 

SIGN 

NORMAL  SURFACE  VELOCITY 
V (REAL)  V„  (IMAG.) 


49  53  57 


data  INPUT  FORM  (2) 


DATE 


FAR  FIELD  PRESSURE  CALCULATION 
TITLE  SAMPLE  PRoBkSn  J 

PROBLEM  


SHEET 


3 


Of  3 


AUTOMATIC  PAR  FIELD  POINT  GENERATION 
NPPLAT^PptNG  LATLIM  LONGLIM 


0/?lQO«flo6^l  I <70.0 


10 


18 


25 


ARBITRARY  FAR  FIELD  POINT  SPECIFICATION 
'^FF  ^F 


*'FF 


XWAVE  DATA 


29 


-*.>0»*E-ll  -J. •••«••• 


9.9‘»52E-0l  ?,IIOJE*QO  J.^90?e-0^  l.a«i54E-0l  9.9»»5ie-«l 


rstSE^Ol  ?.3r9lE-01  -1 .?^Oa£»Oa  2.ai08C*00  S.H907E-02  9.76t9£-01  2.079l£-fit 


I 


39 


43 


» A 


%r2lt-0i  S.O2t2e*01  J.2S9%e«01 


j.i>i42E-gi  a.Jisi-i-oi  -«.<>i92e-oi  •.jjiie-oi  -s.52<i4e-oi  r.2<ii9E-gi 


SuRFAi.  RRESSJR£S<  PART.  IMAGINARY  PART 


50 


51 


33fj£.Q^  9.98‘»3£«ll  -1. 05<*2  i *11 1 8»8bl%E*02  9.bbl5E-0l  9«5li5C*8l 


S.kJ%IE»00  S.SJia&'ll  f.i«7SE-tt  ^.bS^OE^tl  l.lJbSE^Ob  b•J^SlC♦0• 


C3.«riTU0t  > l.SOO£»ll 


L0NCITJ3E  10.  3£.*l  ( ) . 00  B £ *0 1 OEG)  1 1 . 600E«02D£&I  ( 2 • 70  0£ » 020C& I 


1 

L 


58 


f.  % ■ 


PIECEWISE  CONICAL  SHELL  SURFACE 


1 4 

REGION.  EXTENT  AND  MODELING  _ 

'i  '1  *1 

^2 

^2 

n 

m 

SIGN 

1 9 17 

NORMAL  surface  VELOCITY 
V„  (REAL)  V.  (IMAG.) 


1 e 11  18 

REGION:  EXTENT  AND  MODELING  - 

r, 

»2 

Zt 

Z2 

n 

m 

SIGN 

49  53  57 


NORMAL  surface  VELOCITY 
V (REAL)  V„  (IMAG.) 


Fine  iT  E*lT.Fr  IVH  K4£W;  ■ Sf^OMOS 


L. 


61 


63 


FAWICLO  PRESSURES  *T  THE  SURFRCE  OF  A LARGE  SPHERE  CENTERED  AROONO  THE  BOOT 


real  tlAClNARV  REAL  INAGINART  REAL  INAGINART  REAL  IMAGINART 


coLiriTuoE  > i.gagE*oi  dec 


COUTITUOe  = 5.<.0I)E»I)1  DEG 


COL*TITUOC  » OEG 


REAL  r*4AGrNART  REAL  IHAGMART  REAL  IMACIHARt  REAL  XHA6INARY 


LONGITUDE  10.  DEGI 


COLfeTITUOC  a 1.060E»02  DEC 


COLfcTITUOE  s l.?OQE»0?  O^C 


coLATtruof  3 i«}Qor«o2  oec 


86 


COL*TITUD€  » t•(•00E»0^  0€G 


COLATITUOC  3 l.^?0e^0Z  DEG 


COUTITUOE  » l.f>?l)E>92  DEC 


C0L4TITU0E  * l.b«*OE*0?  OEG 


6rt5£-Bt  S.S9ftOF*00 
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1 
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1 
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